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I. INTRODUCTION 

Quantum spin liquids (QSLs) are symmetry-preserving 
quantum phases of matter hosting fractionalized 
excitations.^ They not only possess unusual physical 
properties, which makes them interesting from the per- 
spective of fundamental physics, but may also provide 
new opportunities for quantum information technology.^ 
The S = 1/2 kagome Heisenberg antiferromagnet 
(KHAF) has emerged as a promising candidate system 
from a decades-long search for QSL in theoretical models 
and real materials. Kagome is a lattice of corner-sharing 
triangles. Spins, residing in the vertices of kagome and in- 
teracting antiferromagnetically, are frustrated by the in- 
ability to minimize all pair wise interactions. The strong 
frustration suppresses spontaneous symmetry breaking 
and thus opens the possibility of a spin-liquid state. ^ The 
S = 1/2 KHAF is realized in several materials, most no- 
tably in herbertsmithite ZnCu3(OD)6Cl2, which shows 
signs of a QSL.^ 

The Hamiltonian of the S = 1/2 KHAF is 



{ij) 



(1) 



where are S = 1/2 spin operators defined on vertices 
of kagome, and the summation is over nearest-neighbors 
(Fig. [1^). Despite its simplicity, determining the ground 
state of Eq. ([T]) has proved to be a formidable task. Sev- 
eral competing proposals for the nature of its ground 
state have been made over the past two decades ."^^ 

Recent numerical studies based on the density-matrix 
renormalization group (DMRG) method cast new light 
on this problem. They provide compelling evidence for a 
QSL. In the ground state, spins have very short-ranged 
correlations, gapped S = 1 exci tations, and no signs of 
long-range magnetic order.CSIIIl xhe ground-state wave 
function bears topological entanglement entropy with to- 
tal quantum dimension D = 2, hinting at a spin liquid 
of the Z2 type.l^I^ Additionally, the numerical studies 
have revealed a wealth of new features that remain to be 
understood, including a strong valence-bond resonance 
around diamond-shaped loop^^ and peculiar valence- 



bond correlations near lattice defects^. The ground 
state of Eq. [l] with cylindrical geometry shows a strong 
dependence on the circumference and the chirality of the 
cylinder. Localized S = 1/2 spins are found on the 
open edges of certain types of cylinders .^^ These numer- 
ical findings call for a theoretical assessment. 

On the theory front, much effort has been made to 
construct a low-energy effective description of the Z2 
spin liquid phase of the 5' = 1/2 KHAF. A Z2 spin- 
liquid theory is commonly obtained by mapping a quan- 
tum dimer model (QDM), itself a low-energy description 
of the S = 1/2 Heisenberg model, to a Z2 gauge the- 
ory. This method was initiated by Moessner et al^^ 
and has been successfully applied to the Z2 spin-liquid 
phase of the triangular antiferromagnet^^ and to the Ji- 
J2 Heisenberg model on a square lattice. - - On kagome, an 
exactly solvable QDM with a Z2 spin-liquid ground state 
was constructed by Misguich et al^. Its Z2 fiux (vison) 
excitations are gapped and completely localized. Dynam- 
ical visons were introduced by Nikolic and Senthil^^ and 
by Huh et al^^ and various Z2 spin-liquid phases and 
proximate valence-bond crystals were classified on formal 
grounds. Yet, how the aforementioned numerical results 
fit into a concrete Z2 spin liquid theory remains an open 
question. 

In this paper we provide a phenomenology of the Z2 
spin-liquid phase of the S = 1/2 KHAF, which bridges 
the theory and the numerics. We demonstrate that the 
kagome QDM is naturally described in the language of a 
Z2 gauge theory. We construct a Z2 lattice gauge theory 
and show that the model indeed reproduces qualitatively 
various features of the ground state of the S = 1/2 KHAF 
revealed by the DMRG studies. 

The paper is organized as follows. In Section |TI| we de- 
scribe the construction of the phenomenological model. 
In Section HI we describe techniques for solving the phe- 
nomenological model. In Section [IV| we solve the model 
in various settings and compare the solutions with the 
numerical results. In Section [Vl we discuss the connec- 
tion between the current work with previous works. 
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FIG. 1: (a) The kagome Heisenberg antiferromagnet. S — 1/2 
spins are located on the lattice sites, and the antiferromag- 
netic Heisenberg interactions are between nearest-neighboring 
spins, (b) The rules of the arrow representation of the dimer 
coverings, (c) An example of the arrow representation, (d) 
Mapping the arrow patterns to Ising variables. Shaded (un- 
shaded) triangles are A (B) sites of the honeycomb. 



II. CONSTRUCTION OF THE MODEL 

In this section we construct a phenomenological model 
describing the Z2 spin liquid phase of the S = 1/2 
kagome antiferromagnet. Our starting point is a QDM 
on kagome that is thought to represent low-energy states 
of the S = 1/2 KHAF in the sector with zero total spin. 
We first show that the Hilbert space of kagome QDM 
is identical to that of a Z2 gauge theory on a dual lat- 
tice, a honeycomb whose sites are centers of kagome tri- 
angles. We recast the QDM Hamiltonian of Misguich, 
Serban, and Pasquier (MSP), known to have a Z2 spin- 
liquid ground state, as a Z2 gauge theory. We then per- 
turb the MSP Hamiltonian by adding dimer interactions 
consistent with the lattice symmetry. We keep only the 
leading terms beyond the MSP model. The resulting Z2 
gauge theory is mapped onto a dual Ising model in a 
transverse field, which is used to compute the density 
and correlations of quantum dimers. 



A. From quantum dimer model to Z2 gauge theory 

The starting point is a QDM on kagome introduced 
by Zeng and Elser.!^ A dimer represents two adjacent 
S = 1/2 spins in a state with total 5' = 0. The QDM 
Hamiltonian acts on the space of dimer coverings, states 
in which each site forms a dimer with one nearest neigh- 
bor. The basic assumption of the QDM is that dimer 
coverings span the low-energy subspace of Eqjl] with to- 
tal spin 5* = 0. Elser and Zeng^^ mapped dimer cover- 
ings on kagome to patterns of arrows connecting centers 



of adjacent triangles (and thus forming a honeycomb). 
On a triangle with a dimer, one arrow points out of the 
triangle and away from the dimer; the other two arrows 
point into the triangle. A triangle without a dimer has 
three arrows pointing out. The mapping is illustrated in 
Fig. [T]3. It is easy to check that the arrows are subject 
to a constraint: each vertex of the honeycomb lattice has 
either one outgoing and two incoming arrows, or three 
outgoing arrows. An example of such a mapping is shown 
in Fig. [T]^. (This convention, used by Misguich et al^^ 
is in fact opposite to the one used by Elser and Zeng^^.) 

We introduce Ising variables cf^^ = ctj^ = ±1, defined 
on the links of the honeycomb lattice, to parametrize 
the states of the arrows. To this end, we partition the 
honeycomb sites into A and B sublattices. For a given 
honeycomb link (ij), afj = 1 if the arrow points from 
an A site to a B site, and —1 otherwise (Fig. The 
constraint on arrows is thereby translated to a constraint 
on Ising variables: 



Qi 



1 ieA 
-1 ieB 



(2) 



Here i denotes a site on the honeycomb lattice, whereas 
(il),(i2), and {i3) are the three links emanating from the 
site i. 

We interpret the variables afj as the electric flux op- 
erators of the Z2 gauge theory on a honeycomb. The 
constraint Eq.[2]then becomes Gauss's law, and Qi is the 
Z2 charge on a honeycomb site i. 

We have thus established that the Hilbert space of a 
kagome QDM is identical to that of a Z2 gauge theory on 
a honeycomb with staggered background charges As 
a result, any kagome QDM Hamiltonian can be written 
as a Z2 gauge theory Hamiltonian. 



MSP Hamiltonian as Z2 gauge theory 
Hamiltonian 



The MSP Hamiltonian describes the simplest, exactly 
solvable QDM on kagome.^^ Its ground state is an equal- 
amplitude superposition of all possible dimer coverings, 
the analog of the Kivelson-Rokhsar state^^ on kagome. 
Here we will motivate the MSP Hamiltonian from the 
gauge theory perspective. 

In Section |II A we have introduced the electric flux 
operator afj and related it to the arrow orientation on 
link (ij). Here we consider the Z2 magnetic flux operator 
associated with a hexagonal plaquette of the honeycomb 
(Fig.§i): 



= n 



12^23^34^45^56^61 ' 



(3) 



and erf. are the 



Here 1, 2, ... 6 are sites in hexagon a 
standard Z2 link variables. commutes with charge op- 
erators [F^^Qi] = 0, Va,z. Therefore, the Hilbert space 
subject to the charge constraint Eq. |2]is invariant under 
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FIG. 2: (a) The flux is defined as the product of a^j around 
the hexagon, (b) The magnetic flux threading any hexagon 
is 1 in the ground state, (c) A vison state contains a — 1 
flux threading through a plaquette (shaded hexagon). (d) The 
effect of on the electric fluxes is to reverse all the arrows 
around the hexagonal plaquette a, or equivalently flipping the 
dimers in a David Star. 

F^. Electric flux operators af^- and magnetic flux opera- 
tors F^ are the building blocks of a Z2 gauge theor}^. 

Consider a simple Z2 gauge theory Hamiltonian, con- 
taining only magnetic terms: 

H = -hY,K, (4) 

Oi 

The summation is over all hexagonal plaquettes; h > {). 
F^ reverses all arrows around the hexagon a regardless 
the initial orientation. In terms of dimers on kagome, 
F^ shifts dimers along a closed path within a David star 
with amplitude —h (Fig. [2]i). which is identical to the 
operator ''a^{hy^ defined by Misguich et al.^. Thus, the 
Hamiltonian Eq. |4] is nothing but the MSP Hamiltonian 
formulated in the Z2 gauge theory language. 

The MSP Hamiltonian Eq. [4] is solvable. The eigen- 
states are labeled by the magnetic flux threading through 
every plaquette. The ground state contains no magnetic 
flux, and it has been shown that the ground state re- 
spects all symmetries of the lattice and exhibits short- 
ranged correlation, as expected for a gapped spin liquid 
(Fig. The excited states are created by inserting 

— 1 fluxes, known as visons in literature, and the energy 
cost for each vison is 2h (Fig^). Visons are completely 
localized and dispersionless in the MSP Hamiltonian. 



C. Building the phenomenological Hamiltonian 

The MSP model is unfortunately too simple to account 
for the rich features observed in DMRG numerics. For 
instance, the response to a lattice defect is conflned to 
its immediate neighborhood as opposed to the ripple- 
like features observed in DMRG calculations.^^ Yet, the 




FIG. 3: (a) The links highlighted in red (thick solid) are a pair 
of nearest neighbors, and the links highlighted in blue (thick 
dashed) are another pair, (b) Similar to (a), the links high- 
lighted in the same format are second neighbors, (c) A pair of 
third neighbor links are highlighted in red (thick solid), (d) 
The two quasi-degenerate dimer coverings in the "diamond" 
resonance observed by DMRG calculations, (e) The arrow 
representation for the diamond resonance. 

MSP Hamiltonian provides a natural starting point for 
the construction of a phenomenological Z2 gauge theory. 
In what follows, we perturb the MSP model by adding 
electric flux operators to the Hamiltonian. From the 
quantum dimer model perspective, this amounts to in- 
cluding dimer density interactions. 

We consider generic interactions up to the second order 
in the electric fluxes 

(ij) {ij),{kl) 

The coupling constants ^^ij^i i^iust be compatible 

with the symmetries of the kagome and the honeycomb 
lattices. 

We flrst show that the linear terms are not allowed 
because they are odd under some lattice symmetries. A 
TT rotation about a kagome vertex shared by triangles i 
and j reverses the arrow on honeycomb link (ij). This 
alters the sign of the electric flux, afj hence 

a[^^ = 0. Physically, a linear term expresses a preference 
for some dimer conflguration, any of which violates the 
full symmetry of kagome. 

We next establish that second-order interactions van- 
ish for honeycomb links that share a site, e.g., (12) and 
(24) (nearest neighbors, distance 1/2 lattice spacing) 
shown in Fig. [3^. A second-order term 0-^2^24 can be 
transformed into a flrst-order one, icrfs, by using Gauss's 
law Because linear terms are forbidden, we conclude 

that -^x2'24 ~~ 

By combining Gauss's law with lattice symmetries, we 
can rule out interactions between second neighbors (dis- 
tance V3/2). In Fig. [Id, sites 2 and 5 carry opposite 
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Z2 charges, hence ol^ol^cfl-, = -^^53^56^52. ^21^24 



— (753 crfg. Since = ±1, we obtain crficrfg = — 0-240- 

f2) 

By lattice symmetry, ^21.53 = 
neighbor terms cancel out: 



'24^56- 

(2) 

^24-56 • second- 



j (2) X _x \ aK^^J _x \\^} ( _x _x \^x _x \ fk 

^21;53^21^53^^24;56^24^56 — ^21;53 V^21^53^^24^56>' — 

Thus there are no interactions between second- neighbor 

(2) (2) 
links and we may set ^21-53 — ^24-56 — 0- 

(2) 

The same line of argument shows that ^21-56 ~ 



i(2) 



(2) 



i(2) 
^24;53 



0. This rules out a a^a^ term for links that 
are third neighbors (distance 1) residing on different 
hexagons of the honeycomb. 

The closest non-vanishing interactions are between 
third neighbors (distance 1) residing on the same 
hexagon, e.g., links 24 and 35 in Fig. (Sj). Adding these 
interactions takes us one step beyond the MSP model. 

To determine the sign of the coupling constant between 
these third neighbors, we use input from numerical stud- 
ies. The DMRG calculations of Yan et al^ have iden- 
tified a strong valence-bond resonance in the diamond- 
shaped loop, which corresponds to the tunneling between 
two quasi-degenerate dimer coverings (Fig. [sji). Trans- 
lating the two configurations into electric flux variables 
(Fig. [3^), we see that electric fluxes cr^ on opposite sides 
of a hexagon tend to be opposite, which indicates a neg- 
ative coupling between the electric fluxes. 

We have thus obtained a phenomenological Hamilto- 
nian describing a Z2 spin-liquid state of the kagome 
Heisenberg antiferromagnet: 



H = -hJ2F^ 



X _x 



(6) 



3.n. 



Pairwise interactions of electric fluxes are limited 
to honeycomb links facing each other across a hexagon 
(Fig. [s]^). Physical states satisfy Gauss-law constraints 
K >0, and h > 0. 

III. SOLVING THE MODEL 
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FIG. 4: (a) The dual Ising variables on defined on the hexago- 
nal plaquettes, and they form a triangular lattice (blue dashed 
lines). Latin (Greek) letters are labels of the sites in the hon- 
eycomb (triangular) lattice, (b) The constraints on the phase 
factors Xai3- The product of the Xaf3 around the shaded (un- 
shaded) triangles is —1 (1). A possible arrangements for X^^ 
is that Xa/3 = 1 on black solid links and —1 on red dashed 
links, (c) The four sublattices of the triangular lattice. A, 
C, and D are labels of the sublattices. Sites in the same sub- 
lattice are connected by the third-neighbor interaction, (d) 
The exchange interaction between the spins in A sublattice. 
The interaction is K on black solid links and —K on red 
dashed links. The primitive vectors ^1,2 are also shown. 



translates into a transverse field coupled to new Ising 
variables r^. The interaction of electric fluxes translates 
into a standard Ising coupling r^r| between third neigh- 
bors of the new lattice. The system thus consists of four 
entirely decoupled sublattices (A, B, C, and D in Fig.[4]3) 
and behaves as four decoupled Ising models. 



In this section, we present the techniques needed for 
studying the new model. We first map the Z2 gauge the- 
ory (|6| to a quantum Ising model through a standard 
duality transformation. Details of the duality transfor- 
mation are provided in Section III A In Section III B we 



construct an analytically solvable soft-spin version of the 
dual Ising model. In Section [lII C| we translate Z2 electric 
fluxes in the honeycomb lattice to the density of dimers 
on kagome. The latter is directly related to the nearest- 
neighbor spin correlations in the S = 1/2 kagome Heisen- 
berg antiferromagnet, connecting our phenomenological 
model to the DMRG results. 

The dual Ising model is defined on a new, triangular 
lattice, whose sites are centers of honeycomb plaquettes. 
Fig. [4^. We label the sites of the triangular lattice by 
Greek indices a, /3, 7 ... As usual, the magnetic- flux term 



A. Duality transformation 

The standard procedure to solve a Z2 lattice gauge the- 
ory Hamiltonian is to map it to a quantum Ising model 
through a duality transformation.^^ The electric flux op- 
erator afj is expressed in terms of new Ising variables 
defined on honeycomb plaquettes: 



_X \ n-^ n-^ 



(7) 



where a and j3 are adjacent honeycomb plaquettes shar- 
ing honeycomb link {ik) (Fig. [4^). Aq,/3 = ±1 is a Z2 
phase factor to be defined below. The magnetic flux op- 
erator translates into 



pz ^ ^x^ 



(8) 



Ising operators and r| obey the standard anticommu- 
tation relation {r^^r^} = Sa^- 

The phase factors A are required in the presence of 
nontrivial background charge. Substituting Eq. [7| into 
Gauss's law ([2| expresses the background charge on hon- 
eycomb site a/37 shared by the honeycomb plaquettes a, 
/3, and 7 to a product of A factors: 



a. 



+1 down triangle, 
— 1 up triangle. 



(9) 



Figj4]3 shows a choice of Xc^j^ satisfying Eq. [9j 

Equipped with the duality transformation, we are 
ready to map the Z2 gauge theory Hamiltonian on the 
honeycomb lattice to a quantum Ising model on the dual 
triangular lattice. Substituting Eq. [7| and Eq. [S] into 
Eq. [6) we obtain the dual Hamiltonian, 



^ Ac^^A^/^r^r^. (10) 

;,/3)G3.n. 



The second summation is over third-neighbor dual spin 
pairs. 7 is a triangular-lattice site located between third 
neighbors a and j3. Since the interaction is only among 
third neighbors, the dual triangular lattice splits into four 
non-interacting sublattices A, B, C, and D (Fig.|4j:). The 
Hamiltonian is further simplified. 



H = Ha + Hb + He + Hd • 
The Hamiltonian for the A-sublattice is given by 



(11) 



HA=-hY,T^ + KY, tItI^s^ +kJ2 

reA reA reA 

-i^^ (12) 

reA 

where ^1^2 are primitive vectors generating the A- 
sublattice (Fig. |4]i). Hb,c,d have the same form as Ha- 
Note that Ha is unfrustrated because each triangle con- 
tains two antiferromagnetic links and one ferromagnetic 
link. 

The decomposition of the dual triangular lattice into 
four sublattices is not accidental. It can be shown 
that the four sublattices remain decoupled if further- 
range two-spin interactions are included in the dual Ising 
model. However, the sublattices can be coupled by four- 
spin interactions. This feature of the dual theory is a 
manifestation of the GI/(2, Z3) projective symmetry iden- 
tified by Huh et a/.l^, and we shall come back to this 
point in Section [V| 

The dual Ising model, being unfrustrated, possesses 
two phases at zero temperature. When the dual 



Ising model Eq. 12 is in the paramagnetic phase, which 
corresponds to the Z2 spin liquid phase of the quantum 
dimer model on kagome. When <C i^, it is in the 
magnetic phase, which corresponds to the valenc e bond 
crystal phase of the quantum dimer model.'^^'^ We fo- 
cus on the paramagnetic phase of the dual Ising model 
throughout this paper. 




FIG. 5: (a) An up-triangle on kagome. ctqi, (Jq2i and ctos 
are Z2 electric flux operators on the links of the honeycomb 
lattice. (ii2, <ii3, and d23 are the dimer density modulation 
on the links of kagome. (b) A triangle occupied by a dimer. 
(c) An empty triangle. 



B. Soft spin model 

The dual Ising model Eq. [TT] and Eq. [12] cannot be 
easily solved. We construct an analytically solvable soft 
spin model that is expected to work well deeply in the 
disordered phase of the dual Ising model."^ To this end, 
we replace the dual Ising operator = ±1 by a real- 
valued operator (pr- The interaction between dual Ising 
spins therefore becomes 

reA 

The transverse field term —hr^ brings quantum fiuctua- 
tions to the dual Ising spin, so we add 



(13b) 



reA 



reA 



to the Hamiltonian. tt-t is the canonical momentum op- 
erator conjugate to 0^. Here a mass term with A > is 
introduced to ensure the fiuctuations of (j)r are gapped. 

Gombining all the terms listed above, we obtain the 
soft spin Hamiltonian for the A-sublattice, 



reA reA 



reA 



K ^ (/)r(/>r+52 ~ ^ (/>r(/>r+6i-h62 (14) 



reA 



reA 



The soft spin Hamiltonian for sublattices 5, C, and D 
are obtained in the same way. The total soft spin Hamil- 
tonian is 



H 



soft 



H 



soft 



H 



soft 



H' 



soft 
D • 



(15) 



C. Dimer density 



The last piece of our construction is a formula relating 
the expectation value of electric fiux operators to the 
density of dimers on kagome. To start with, we consider 
an up triangle on kagome, which corresponds to an A 
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site in the honeycomb lattice (Fig. |5p). The case for a 
down triangle or a B site is obtained by the same token. 
The dimer density modulation on kagome link (ij), dij^ 
is defined as the dimer density on the said link minus the 
average number of dimers per link, 1/4. 

We postulate the following relation between the dimer 
density modulation and the expectation value of Z2 elec- 
tric fiux operators on symmetry ground (Fig. [sfe 



(16) 



where li, v are coefficients to be determined. The ex- 
pressions for di2 and ^13 are obtained by 120° rotations 
around the center of the triangle. 

We consider the case in which one side of the triangle 
is occupied by a dimer and the corresponding arrow con- 
figuration (Fig.jsjD). On the one hand, the dimer density 
modulations are given by. 



112 



lis 



^23 



(17) 



On the other hand, the expectation values of the Z2 elec- 
tric fiux operators are given by 



-1; 



«3) 



-1. 



(18) 



The postulated relation Eq. 16 is consistent with the 
above if and only if. 



1 1 

u=-; v = --. 

Therefore, we obtain the following relation: 
1 



23 



(('^01) - ('^02) - ('^os))- 



(19) 



(20) 



The self-consistency of Eq. 20 can be further examined 
by considering another configuration shown in Fig. (Sj), 
where di2 = c^i3 = c^23 = "i, and (a^^) = {(7^2) 
((703 ) = 1- It can be seen that Eq. 20 holds in this special 
case as well. 



IV. RESULTS 



In Sections O and III 



we have constructed a phe- 
nomenological Z2 lattice gauge theory and developed nec- 
essary techniques to solve it deeply in the deconfined 
phase. In this section we compare predictions of our 
model with the DMRG results in various settings. 

In Sections [lVAl|IVB[ and |IVC[ we study the ground 
state of the phenomenological model on kagome cylin- 
ders. We show that the ground state strongly depends on 



the cylindrical geometry. In Section [IV D[ we study the 
singlet-pinning effect in the ground state of kagome cylin- 
ders. In Section pV E[ we analyze the effect of the open 
boundary on the ground state of kagome cylinders and 
show that certain types of open boundary bind spinous 
near the edge. In all of the cases studied, we find good 
agreement between the phenomenological model and the 
DMRG. 



A. Ground state of YC4m kagome cylinders 

Most of the recent DMRG studies on 6* = 1/2 kagome 
Heisenberg a ntiferro magnet are performed with cylindri- 
cal geometry.flS^im^ The periodic boundary condition is 
imposed in one direction of kagome and open boundary 
condition in the other direction, making the system ef- 
fectively a cylinder. The circumference of the cylinder 
is usually much smaller than the length, and we shall 
assume the length is infinite in the following discussion. 

Th e DM RG study has identified an even-odd 
effect .flS^ A kagome cylinder can be regarded as a one- 
dimensional system, and the unit cell is defined as the 
smallest building block that generates the whole cylin- 
der by translation in the length direction. The ground 
state is a uniform spin liquid when the number of spins 
per unit cell is even. When the number is odd, the spin- 
liquid ground state coexists with a valence-bond density 
wave pattern, which breaks the spatial symmetries of the 
cylinder. 

Motivated by the DMRG findings, we study the phe- 
nomenological model on cylinders. To start with, we con- 
sider the so-called "YC4m" cylinders, where m G Z and 
the number 4m denotes the circumference of the tube.^ 
Note that the number of spins per unit cell is 6m, which 
is even. Fig. |6^ shows a YC8 cylinder. We shall use the 
YC8 cylinder as an illustrative example in the following 
discussion, and all the results can be easily generalized 
to all YC4m cylinders. 

The corresponding Z2 gauge theory model is defined 
on the honeycomb lattice with "armchair" type periodic 
boundary condition (Fig. . Note that the honeycomb 
cylinder is unfolded in a way different from Fig. [6^. 

We introduce a labeling system that shall be employed 
in later discussions. The honeycomb plaquettes are gen- 
erated by primitive vectors e^^ and e^. The position vec- 
tor of a honeycomb plaquette (or a dual triangular site) 
is uniquely expressed as u^u + ve^^ and it is labeled as 
(ix, v). Then, the three honeycomb bonds adjacent to the 
said plaquette are labeled as {u^ a), {u^ 6), and (u,v,c) 
(Fig#). 

The Z2 gauge theory Hamiltonian possesses a 
global conservation law associated with the cylindri- 
cal topology.^ We consider a non-contractible contour 
around the cylinder, shown as thick red bonds in Fig.|6]:3. 
We define the total electric fiux passing the contour. 



X 



(21) 



It can be shown that X commutes with the Hamiltonian 
[X, H] = 0, meaning the total Z2 electric fiux is con- 
served. Therefore, the Hilbert space is divided into two 
topological sectors, with X = 1 and X = —1 respectively. 

The dual Ising model is defined on triangular lattice 
with cylindrical geometry (Figs. |6]i,e). The global con- 
servation law on total Z2 electric fiux X gives rise to 
new constraints on the phase factors Xa(3 in addition to 
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(a) 12 3 4 




(f) ii = 0u = 2u = 4 (g) u = 0u = 2u = 4 

V = 4 A A A V = 4 A A A 

V = 2 A — A — A ^ V = 2 A — A — A — ^ 

I \ I \ I \ I \ I \. I \ 

V = A — A — A ^ V = A — A — A ^ 

FIG. 6: (a) The YC8 cylinder. The periodic boundary con- 
dition is imposed along the vertical direction. Sites with the 
same number are identical, (b) The labeling system employed 
in the main text, (c) The corresponding Z2 gauge theory 
model on the honeycomb lattice. Arrows are the primitive 
vectors eu and of the dual triangular lattices. The coor- 
dinates of the dual triangular sites are shown, (d) The dual 
model for YC8 cylinder in the even sector. The phase factor 
A = 1( — 1) on black solid (red dashed) links, (e) The same 
dual model in the odd sector, (f) The exchange interaction 
between dual spins in A sublattice in the even sector. The 
interaction is K on black solid links and —K on red dashed 
links, (g) The exchange interaction between dual spins in A 
sublattice in the odd sector. 



Eq. [9| Substituting Eq. [7| into Eq. [21] we rewrite the to- 
tal Z2 electric flux operator X in terms of dual variables 
(Fig|6}3): 



(1,^ + 1) 



±1 



(22) 



Thus, the total electric flux X is dual to the total flux of 
Xai3 around the cylinder circumference. 

In what follows, we treat the even {X = 1) and odd 
{X = —1) sectors separately. We consider the even sector 



at first. A choice of Xaf3 satisfying conditions Eq. loland 
Eq.[22]is shown in Fig.jGj). We see that the dual model is 
decomposed into four independent copies of unfrustrated 
quantum Ising model, corresponding to four sublattices 
of the triangular lattice. In the disordered phase of the 
dual model, or equivalently the deconfined phase of the 
Z2 gauge theory model. 



K-) = Xa^ir^rl) = Xap{r^){Tl) = 0. 



(23) 



The second equality follows from the fact that dual spins 
a and /3, being nearest neighbors to each other, must 
belong to two different sublattices. The third equation 
follows from the assumption that the dual model is in the 
disordered phase. Given the modulation of dimer density 
is proportional to (cr^), we deduce that the ground state 
is uniform for YG4m cylinders in the even sector. 

We proceed to calculate the ground state energy of the 
YC4m cylinder in the even sector. The circumference 
of each sublattice is m. The soft spin model for the A 
sublattice is given by (Fig. |6]F) 



rjsoft 



(24) 



(n,'t;)GA (n,i;)GA 



Here the periodic boundary condition (j)u^v = ^u,v-\-2m is 
imposed. The above Hamiltonian can be readily diago- 
nalized. The ground state energy per dual site is given 

by 



m-1 .27r 



dku /, 27rv\ 
27r \ m J 



(25) 



where 



= Y^A + 2i^[cos ku + cos ky — cos(/c^ — ky)] (26) 



is the dispersion relation of soft spin fluctuations. 

The analysis on the dual Ising model in the odd sector 
can be carried out in the similar manner. A choice of Xai3 
is shown in Fig. [6^. Similar to the even sector, the dual 
triangular lattice is decomposed into four independent, 
identical sublattices, and the ground state is uniform in 
the odd sector as well. 

As shown in Fig. |6^, the explicit translational invari- 
ance of the sublattice Hamiltonian is lost in the circum- 
ference direction. However, we can make a gauge trans- 
formation 0n,2m — ^n,2m to rcstorc the cxpHcit trans- 
lational invariance. The soft spin model becomes iden- 
tical to Eq. [24j The boundary condition, however, be- 
comes anti-periodic: (j)u,v = — 0n,v+2m- The grouud statc 
energy per site is given by 



m-l .27r 



e-1 



dk^ ( (2^ + l)^ \ 

-— ku, , (27) 

27r V m ' 



Here uj is the same as Eq. 26 
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FIG. 7: Ground state energy per site e of YC4m cylinder in 
even X = 1 (blue crosses) and odd X = — 1 sectors (red open 
circles). The solid and dashed lines correspond to K/ A = 0.12 
and K/A = 0.15, respectively. 



FigjT] shows the ground state energy per dual site as 
a function of cylinder circumference for generic values of 



A and based on numerical evaluation of Eq. 25 and 



Eq. 27 We find that the ground state energy in the two 
sectors becomes degenerate in the limit of infinite circum- 
ference as expected for a gapped Z2 spin liquid. When 
the circumference is finite, the energy splitting shows an 
alternating pattern: the even ground state has lower en- 
ergy when m is even, and the odd ground state has lower 
energy otherwise. Such a pattern plays an important role 



in our discussion on edge spinous in Section IV E 



To summarize, we have found that the ground state of 
the YC4m cylinder is the uniform Z2 spin liquid. The 
ground state is in the even sector X = 1 for even m, and 
the odd sector X = —1 for odd m. 



B. Ground state of YC4m + 2 kagome cylinders 

In this section, we discuss another family of kagome 
cylinders known as YC4m + 2, m G Z.^^ The number of 
spins per unit cell is 6m + 3, which is odd. Based on the 
Lieb-Schultz-Mattis theorem, we expect that the ground 
state is either gapless or symmetry-breaking.^^ Indeed, 
DMRG has identified a symmetry-breaking valence-bond 
density wave pattern. We will show that our model re- 
produces the same pattern in the ground state. 

As a concrete example, we consider the YC6 cylin- 
der, a member of the YC4m + 2 family, in the following 
discussion. The analysis is carried on in parallel with 
Section [TVAl 



Fig. [8p shows the YC6 kagome cylinder, 

and Fig. [Spshows its corresponding honeycomb cylinder 
on which the Z2 gauge theory is defined. The total Z2 
electric flux in the Z2 gauge theory model is defined as 




(3,2) 
(3,1) 
(3,0) 



v=3 A A A 

I \ I \ I \ 

v = 2 A — A — A ^ 

I \ I \ I \ 

B ^=1 A — A — A 

I \ I \ I \ 

17 = A — A — A ^ 

u = u = 2u = 4 



FIG. 8: (a) A YC6 cylinder and the predicted valence-bond 
density modulation pattern. Lattice sites with the same nu- 
meral label are identical, (b) The corresponding Z2 gauge 
theory, (c) The dual Ising model. A and B are sublattice 
labels. A = 1(— 1) on black solid (red dashed) links, (d) The 
interaction between A dual spins. The interaction is K(—K) 
on black solid (red dashed) links. 



bonds in Fig. |8^. [X, i^] = 0, and the Hilbert space falls 
into two topological sectors, X = 1 (even) and X = —1 
(odd). 

The two topological sectors are related by translational 
operation. To see this, we consider another contour 
shown as the thick green bonds in Fig. ^jp. The asso- 
ciated total electric flux operator is defined as: 



(29) 



On the one hand, it can be seen that X' = —X on ac- 
count of Gauss's law On the other hand, the new 
contour is related to the old one by a e^^ shift, and we 
have T{eu)XT{euy = X' . Here T{eu) denotes the asso- 
ciated shift operator. Therefore, 



T(eOX = -XT(e,). 



(30) 



(28) 



Eq. 30 implies that the two topological sectors are de- 
generate in energy in the YC6 cylinder. Note T(e^) is a 
symmetry. 



corresponding to the contour highlighted as the thick red 



[T(e,), H]={). 



(31) 



9 



Let be a common eigenstate of the Hamiltonian H 
and the total electric flux X\ 



H\iP)=e\i^), X\i,)=x[ 



(32) 



Then, 



HT{eM = T{eu)H\iP) = er(e„)|t/-), 
XTie^m = -r(e„)X|^) = -xT{e^m. (33) 

We see |?/^) and T{eu)\i^) ^^'^ degenerate in energy yet 
belong to two different topological sectors. Therefore, 
the energy spectra in the two topological sectors coin- 
cide, and the degenerate eigenstates are related by the 
translation. Most importantly, in contrast to the YC4m 
cylinders, the ground states in the two topological sectors 
of the YC4m + 2 cylinders are exactly degenerate, and 
they break the translation symmetry in e^^ direction. 

We proceed to discuss the dual soft spin model. We 
shall consider the even sector only (Fig.jsJ:). The circum- 
ference of the dual triangular lattice is 2m + 1, making it 
impossible to partition the lattice into four independent 
sublattices. Consequently, there are only two indepen- 
dent sublattices, labeled as A and B. The Hamiltonians 
for both sublattices are identical. The soft spin model 
for the A sublattice is given by (Fig. P 



rrsoft 



K 



(34) 



The periodic boundary condition ^^^o = ^n,2m+i is im- 
posed. 

We have deduced on symmetry ground that the ground 
state in the even sector must break the translational sym- 
metry. In what follows, we show that the translation 
symmetry is broken by a dimer-density modulation pat- 
tern. To this end, we calculate the expectation value 
of Z2 electric flux in the ground state. For the {u^v^b) 
bonds, 



(35) 



The second equality follows from the fact that the dual 
sites (u^v) and {u -\- l^v) belong to different sublattices 
and therefore are uncorrelated (Fig.jsj^). The third equal- 
ity follows from the assumption that the dual model is in 
the disordered phase. By the same token. 



(36) 



The expectation value of ^ ^ does not vanish in the 
ground state. 



-\-d if u = 0mod2, 
—d if 1^ = 1 mod 2. 



(37) 
(38) 



Xi 10"' 

^10' 
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FIG. 9: Dimer density modulation |d| in the ground states of 
Y 04:771 + 2 cylinders as a function of the cylinder index 771. 



Here d ^ in general because the dual sites (u^v) and 
{u^v + 1) belong to the same sublattice (Fig. [s]^). The 
plus and minus signs come from the A phase factor. We 
have defined 



d = 



(39) 



where d doesn't depend on u^v thanks to the transla- 
tional invariance of the sublattice Hamiltonian and the 
fact that the Hamiltonians for A and B sublattices are 
identical. 

and 



Combing Eq. 35 



38, we conclude 



Eq. |36) and Eq. 

that the Z2 lattice gauge theory displays an alternating 
electric-flux pattern in the ground state on the YC4m + 2 
cylinders. The pattern in the even sector is shown in 
Fig. [Sj: as arrows. The pattern in the odd sector is ob- 
tained by a shift, or equivalently by reversing all the 
arrows. 

The Z2 electric flux operator is tied to the density mod- 
ulation. We can therefore determine the dimer density 
modulation in the ground state of YC4m + 2 cylinders by 
using the results in Section pll C| The dimer modulation 
pattern is shown in Fig. [S^i. (The pattern in the other 
sector is obtained by a shift of e^.) The blue solid and 
red dashed bonds stand for dimer density modulation of 
d/4: and —d/A respectively. In terms of the Heisenberg 
model, they correspond to stronger and weaker nearest- 
neighbor spin-spin correlation. The pattern shown in 
Fig. ^ agrees well with the DMRG result .2^1 

The value of d can be explicitly calculated by using the 



soft spin model. We obtain, 

dku cos[27rv/(2m + 1)] 



d- 



1 



E 



2(2m + l)— ^0 



27r uj[ku,2'Kv/{2m + l)] 



, (40) 



where 



cos ku -\- cos 2ky — cos{ku -\- 2ky)]. (41) 



Fig. [9] shows the absolute value of d as a function of 
cylinder index m for generic value of model parameters. 
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FIG. 10: (a) Valence-bond modulation pattern on the YC9-2 
cylinder calculated from the phenomenological theory. The 
spin-spin correlation is enhanced on the solid blue bonds and 
weakened on the dashed red bonds, respectively. Sites labeled 
with the same number are identical, (b) The corresponding 
gauge theory model and the dual Ising model. Arrows indicate 
the non- vanishing expectation value of afj in the ground state. 
Boundary bonds labeled with the same number are identical. 
Thick green bonds highlight the contour with respect to which 
the total electric flux operator X is defined. A and B stand 
for two independent sublattice sites of the dual lattice. The 
phase factor A = 1( — 1) on black solid (red dashed) bonds. 
Here the even sector X = 1 is shown. 



We can see that \d\ decreases almost exponentially as m 
increases, reflecting the finite correlation length in the 
ground state. 

To sum up, we have found that the ground states of 
the YC4m + 2 cylinders are two- fold degenerate. The 
ground states display a symmetry-breaking valence-bond 
modulation pattern, which agrees with what found by 
DMRG. 



C. Ground states of YC(4m+ l)-2 kagome cylinders 

In this section, we briefly discuss one more family of 
kagome cylinders dubbed YC(4m+l)-2^^. This new fam- 
ily of cylinders lacks for the reflection symmetry along 
the circumference direction. The DMRG calculation has 
identified an valence-bond modulation pattern similar to 
the one found in the YC4m+2 family but with a different 
orientation. The analysis on YC(4m + l)-2 cylinders can 
be carried out in the same manner as in Section IV B[ 
and therefore we simply summarize the result. 

Fig. [TO^ shows the YC9-2 cyliner, a member of the 
YC(4m + l)-2 family. Note that the periodic bound- 
ary condition along the circumference direction of cylin- 
der is accompanied by a shift in the length direction. 




FIG. 11: (a) A singlet or dimer is pinned on a bond with en- 
hanced Heisenberg exchange interaction, (b) A pinned dimer 
is amount to frozen arrows in the arrow representation, (c) 
The dimer density modulation pattern on YC8 cylinder cal- 
culated from the phenomenological theory. We choose A = 1, 
K = 0.15, and V = 0.2. The dimer density is increased on 
blue solid bonds and decreased on red dashed bonds. The 
thickness is proportional to the magnitude of modulation. 



Fig. [TOp shows the corresponding gauge theory model 
and the dual Ising model. The topological sectors are 
defined with respect to the contour highlighted as the 
thick green bonds. Similar to the YC4m + 2 cylinders, 
the energy spectra in both sectors are exactly identical. 
In particular, the ground states are two- fold degenerate, 
and they break the translation symmetry along the 
direction. Here the expectation value of afj in the even- 
sector ground state is shown as arrows, from which we de- 
termine the valence-bond modulation pattern (Fig.jlO^). 
The pattern in the other ground state is obtained by a 
translation. 

Comparing the calculated pattern with the one ob- 
tained by DMRG, we find that the two agree quali- 
tatively. However, the pattern predicted by our phe- 
nomenological theory is perfectly uniform along the e^^ 
direction, whereas the DMRG pattern shows another, 
weaker, modulation along the directiorP^. This dis- 
crepancy may be due to the shorter-range physics that is 
not included in our theory. 



D. Singlet-pinning effects in the ground state 

In this section, we study the response of the Z2 spin liq- 
uid ground state to an externally-pinned singlet. Specif- 
ically, when the Heisenberg exchange interaction on a 
bond is enhanced by amount of 6 J, it is energetically fa- 
vorable to form a spin singlet on that bond (Fig. pT^i). 
The pinned singlet freezes the resonating singlets nearby. 
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and a modulation pattern of nearest-neighbor spin-spin 
correlation (S^ • Sj) appears in its neighborhood. When 
<C 1, the pattern coincides with the dimer-dimer cor- 
relation in the unperturbed ground state. The goal of this 
section is to compute such a modulation pattern by using 
the phenomenological model and to compare it with the 
numerical results. 

The effects of an enhanced bond in the Heisenberg 
model are effectively described by an dimer-attracting 
potential in the quantum dimer model, which is amount 
to introducing favored arrow orientations in the arrow 
representation (Fig|TT]3). The Z2 gauge theory Hamilto- 
nian is given by 



H = Ho- F«o,a - 



0,0,6 



^0,1, o'' 



(42) 



where Hq is the unperturbed Hamiltonian and V is the 
strength of the pinning potential. The above Hamilto- 
nian is readily solved by using the duality mapping and 
the soft spin approximation. 

Fig. \n\ : shows the calculated dimer density modula- 
tion pattern on the YC8 cylinder by solving Eq. 42 for 
A = 1, K = 0.15, and V = 0.2. Note that the unper- 
turbed ground state (V = 0) of the YC8 cylinder is uni- 
form as discussed in SectionjiVAj The perturbed Hamil- 
tonian is solved in the even sector, where the ground 
state of the YC8 cylinder lies (Figj7|. The dimer density 
is increased on the blue solid bonds and decreased on the 
red dashed bonds. The bond thickness is proportional to 
the magnitude of modulation. The bond with enhanced 
exchange interaction has the largest dimer density. 

We find qualitative agreement between the phe- 
nomenological theory and the DMRG.*^ In addition, we 
note the following features of the modulation pattern cal- 
culated from the phenomenological theory. Firstly, the 
spatial modulation quickly decays as the distance to the 
enhanced bond increases, which is another manifestation 
of the gapped nature of the ground state. Secondly, there 
are alternatively enhanced and reduced diamond shapes 
which tile the whole cylinder, highlighted as plus and mi- 
nus signs in FigfTT]^. These diamond shapes coincide with 
the so-called "diamond resonance" observed by DMRG 
and Lanczos diagonalization studies^^. 



E. Edge spinons and Z2 screening 

In our previous discussion, we assume that the kagome 
cylinders are infinite in the length direction. As we shall 
see below, the open boundaries give rise to a new phe- 
nomenon: in certain circumstances, an open boundary 
can bind a spinon nearby. The edge spinons are de- 
tectable in numerical simulation, which could serve as 
a direct evidence for fractionalization.^^ 

To start with, we consider a particular type of open 
boundary of YC4m cylinders, which is employed by Yan 
et alJ^ (Fig- 12 i). The motivation is to make the open 
boundary compatible with the "diamond resonance" pat- 
tern. It is energetically favorable to form singlets on dan- 







X = 1 



X = -1 



1 



FIG. 12: (a) The open boundary of YC4m cylinders employed 
by Yan et al . The singlets are pinned at the boundary, 
(b) The arrow representation for the frozen dimer covering 
shown in (a). The numbers ±1 are the Z2 electric fluxes on 
the dangling bonds, (c) Another open boundary condition 
of YC4m cylinders and the pinned singlets, (d) The arrow 
representation and the Z2 electric fluxes corresponding to the 
frozen dimer covering shown in (c). (e) The electric flux in 
the bulk is screened by a spinon. 



gling bonds, and therefore a particular dimer covering 
is pinned at the boundary, which has been observed in 
DMRG calculations. 

The physical meaning of such an open boundary con- 
dition becomes transparent in the Z2 gauge theory de- 
scription. Converting the frozen dimer covering to arrow 
representation, we find that the effect of the open bound- 
ary is to enforce an alternating Z2 electric flux pattern at 
the edge (Fig. 12 3). The total edge electric flu x is (- 1)^ 
in the YC4m cylinder. As discussed in Section pVB[ the 
total electric flux in the ground state of YC4m cylinder 
is ( — 1)^. Comparing the total electric flux enforced by 
the open boundary condition with the one selected by 
the bulk energetics, we find that the two match for any 
m. 

However, the agreement between the edge total electric 
flux and the bulk total electric flux can be lost if one 
chooses other types of open boundary. An example is 



shown in Fig. 12z. Similarly, dimers are pinned at the 
boundary, and an electric flux pattern is enforced at the 
edge (Fig. 121). The total electric flux at the edge is 
always 1 for any YC4m cylinder. When m is odd, the 
edge total electric flux and the bulk total electric flux 
don't match. 

The mismatching is resolved by a spinon screening the 
bulk electric flux in the ground state of YC4m cylinder 
when m is odd, which is in close analogy to the screening 
phenomenon in a metal (Fig.[l2^). Close to the edge, the 
total electric flux is 1, and the cylinder segment between 
the edge and the spinon is in the even sector. Deep into 
the bulk, the total electric flux is —1, and the bulk is in 
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the odd sector. Assuming that the distance between the 
edge and the spinon is x, the total energy of the system 
is approximately given by: 

£^ ^ [ei(m) — e_i(m)]x + const, (modd). (43) 

Here ei(m) and e_i(m) are energy per unit length in 
the even and odd sectors of the YC4m cylinder. Note 
that ei > e_i when m is odd. Therefore, the spinon is 
confined to the edge by a linear potential. 

Finally, we remark that, as the two topological sectors 
in YC4m + 2 cylinders are exactly degenerate, there are 
no edge spinous in generic conditions. 



V. DISCUSSION 

We close our paper by placing our work in the context 
of previous theoretical efforts. As we have shown in Sec- 
tion II A, the hard core constraints of dimers naturally 
endow the kagome QDM a Z2 gauge theory structure, 
by which any kagome QDM Hamiltonian can be formu- 
lated as a Z2 lattice gauge theory Hamiltonian. Earlier 
workers^^ used an alternative Z2 gauge-theoretic de- 
scription of the kagome QDM. In that framework, the 
dimer occupation number on each kagome link is identi- 
cal to the Z2 electric flux on that link. The dimer hard- 
core constraints are approximately enforced by Gauss's 
law and the energetics. The effective Z2 gauge theory 
is therefore a Z2 gauge theory on kagome with charge 
Q = — 1 on every site. 

Despite the superficial differences, the two frameworks 
are quite similar. In the previously used framework, vi- 
sons are defined on both triangular and hexagonal pla- 
quettes of kagome. A vison is associated with quantum 



tunneling of dimer coverings in a plaquette. However, 
dimer coverings on a triangular plaquette cannot tunnel 
as they are eigenstates of the local interactions, and thus 
it is natural to treat the visons on triangular plaquettes 
as high energy degrees of freedom and integrate them 
out. The resulting theory should be the same as ours. 
Our approach thus provides a more economical represen- 
tation of the QDM on kagome. 

The similarity between the two frameworks becomes 
more evident in the dual formulation. The dual Landau- 
Ginzburg order parameter of Huh et al.^ has four in- 
dependent real components <pi. The only quadratic term 
allowed by the GL(2,Z3) symmetry of the dual theory 
is , which is related to the fact that the four sub- 

lattices cannot be coupled at the quadratic level as we 
discussed in Section ttll A[ 

After submission of this paper we learned of a similar 
work by Ju and Balents,'^who mapped a quantum dimer 
model on kagome onto a Z2 gauge theory on a dice lat- 
tice. They computed the dimer modulations in cylinders 
with various boundary conditions and obtained results in 
agreement with ours. 
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